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Abstract — The geoid is used as a transformation linkage between ellipsoidal heights (h) determined from DGPS 
observations and orthometric heights (H). Widespread acceptability and adoption of GPS in local geospatial 
data acquisitions require the development of a local geoid model (N) for use to obtain orthometric heights in the 
absence of a national geoid model. Geoid model can be developed by gravimetric approach; global geopotential 
model (GGM); geometric technique among others. The conventional approach to GPS measurements is the use 
of one base reference station for field measurements. It has several drawbacks e.g. in signal range/coverage, 
accuracy degradation of results, etc. Based on Grashof s law of stability of triangles, this study was therefore 
based on dual reference base stations to improve on DGPS signal range and stability of results. Pro-online 
matrix solver was applied to the least squares observation equations of the two modelled FCT surfaces (multi - 
quadratic and bicubic) to determine polynomial coefficients. The geoid undulation was computed and 
orthometric height generated for production of a topographical plan at lm contour interval for elevation data in 
surveying, engineering and environmental applications. Skill =1 and bias = 0 were computed to confirm the 
predictive capability of the models and that no bias/errors were introduced into the respective modelling 
exercise. Diagnostic test also confirmed the viability and feasibility of providing vertical datum surface for FCT 
by this approach. Standard deviation (a) as accuracy indicator was computed and the multi-quadratic model 
with a =llcm was the better geoid surface for modelling of orthometric height in the FCT by the geometric 
method. 

Keywords — Geoid undulation, Multiquadratic, Bicubic, Grashof’s law, Orthometric height. 


I. INTRODUCTION 

The use of GNSS in orthometric height (H) 
determination requires a geoid model (N) to transform the 
observed ellipsoidal height (h). For global applications, 
global geoid models (EGM2008) have been developed to 
provide the geoid undulation. For small to medium-sized 
areas, global geoid model, according to Odera and 
Fukuda (2015) is too generalized and will lead to error in 
orthometric height if applied. Merry (2009) gave a value 
of 3 m in Central Mozambique when compared with 
EGM2008 values due to the use of generated gravity 
anomalies. Hence, this requires the development of local 
geoid models for the needs of GNSS user community in 
geospatial data acquisitions and applications. 

Al-kragy et al. (2015) observed that a geoid model is 
a three dimensional (3D) geospatial model that defines 
the relationship between the ellipsoid and the geoid 
surfaces at a specific area. Eteje et al. (2018) defined 
geoid as the surface which coincides with that surface to 
which the oceans will conform over the entire earth if free 
www.iieab.com 


to adjust to the combined effects of earth's mass attraction 
and the centrifugal force of earth's rotation. Methods of 
geoid undulation determination are namely: 

(a). Gravity measurements for gravimetric geoid by 
solving general Stoke's integral formula by spherical 
harmonic expansion as given by Heiskanen and Moritz 
(1967): 

A7 „KSm SW R rr w 

N = K TT~ + ^ AsS(,p) ‘ la m 

S(tp) = 1h ---6 sin —-5cos tp- 

. (p 2 

sin — 

2 (2) 

3 cos cp ln(sin ^ + sin 2 

where the various parameters are as given in the 
literature. 

Assuming that the mass of real earth is equal to the 
mass of the normal earth and the potential generated by 
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two masses to be equal, the first two terms in equation (1) 
become zero i.e. 

N = jj A gS ((p)d<y (3) 

The difficulty with Stake's formula is that solution 
requires gravity data all over the earth which is 
impossible to achieve arising from the double integral in 
the formula. To overcome this, global geopotential 
models (GGM) were developed. These global models are 
inadequate for local applications and Odera et al. (2015) 
stated that they are too generalized to be useful for local 
applications and hence, for areas of limited sizes, a local 
geometric geoid model could be developed for 
orthometric data acquisitions. 

(b) Geometric Geoid Model 

This is developed for areas ranging from small to 
medium and computed directly from GPS based 
ellipsoidal height (h) and collocated with points of known 
orthometric heights (H). From the relationship given by 
Kotsakis and Sideris (1999), Jekeli (2006), a linear 
relationship between h, H and N where E is deviation of 
the vertical and curvature of plumb line is 
h = N + H+% (4) 

Seker and Yildirin (2002) observed that at E, = 1", the 
error incurred is 0.08mm which is negligible, 
insignificant and of no practical consequence. Also, 
Nordin (2009) computed the effect of E, = 1" as less than 
lmm. Figure 1 shows the linear relationship between the 
heights. The combined interpretation and implication of 
the above values is that we can write with confidence 
that: 

N = h-H (5) 



Fig. 1: Relationship between Orthometric, Geoid and 
Ellipsoidal Heights 

Source: Olio (2013) and Eteje et al. (2018) 


1.2 Datum Ambiguity/Bias 

The vertical datum may be inconsistent due to 
measurement error and is evident in all three height 
components of h, H and N. Let 3h, 3N and 3H be the 
biases in each of the components, we can, therefore, write 
from equation (6) 
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h-N + H 

(6) 

that 


h + Sh = (H + SH) + (N + SN) 

(7) 

h = H + N + (SH + SN-Sh) 

(8) 

h- H + N + S 

(9) 

where, 



S = (SH + 5N - Sh) 


By comparing h = N + H + E, with h = H + N + S , 
it can be shown that S = E,. Hence 
i E, = (SH + SN — Sit) is insignificant. From the various 
values of E, computed by the above authors, the datum 
bias can therefore be taken as insignificant and hence 
negligible for low order survey and engineering 
applications and adequate for geometric geoid modelling 
(from N = h — H) and hence orthometric height 
determination from H = h — N. Milbert and Smith 
(1996) observed that the very small values of S compared 
to N support the direct conversion between ellipsoid and 
orthometric vertical datum even if they are not defined on 
a common reference. Geometric geoid model hence is 
adopted for modelling orthometric height in the provision 
of vertical datum for elevation data acquisition. 

Kamaludin et al. (2005) observed that differential 
heighting method can be used to eliminate datum 
inconsistencies for topographical and 

engineering/environmental studies and applications. 
From N —h — H, interpolation of geoidal heights ( N ) 
becomes feasible over interested points with an available 
GPS ellipsoidal and existing orthometric heights. 

1.3 Justification of Adopted Field Procedure 

Generally, in DGPS campaigns only one base 
reference station is adopted for observations in the 
relative approach. This method has limitations in 
coverage and accuracy is spatially degraded after certain 
distances beyond, for example, 10km or over large areas. 
Martensson (2002) recommended the use of network that 
resembles a triangulation network in GPS campaigns 
where the aim is to obtain surface cover for geometric 
geoid modelling to ensure that no deterioration of results 
are experienced and hence it can be stated that the results 
from this study are highly stable and consistent since the 
FCT triangulation network was used for all 
measurements. Chang and Lin (1999) reported from 
studies using one and multiple base reference stations, 
that results obtained from the latter are more reliable and 
consistent achieving over 60% improvement in values 
both in horizontal and vertical components using DGPS. 
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1.4 Stability of Shapes 

A triangle is the simplest of closed figures in two 
dimensions and described as the strongest geometrical 
shape and most stable too because of its inherent 
structural characteristics. For example, a square is capable 
of becoming a parallelogram whereas a triangle is only 
capable of being a triangle. The explanation given to why 
the triangle is more stable than other shapes is that it only 
takes three points to define a plane. By adding any point 
to the plane will make it harder and harder for it to be 
stable. Also, no matter where the vertices of the triangle 
lie, they will always define a plane and hence triangles are 
both stable and rigid. Grashof’s relationship can be used 
to compute geometrical stability of figure from Quora 
http:// www.quora.com as: 

n = 3(L — 1) — 27 — h (10) 

where, 

n = degree of freedom 

L is number of links; J is number of joints ; h is number 
of higher points. If n = 0, there is geometrical stability 
of results for a triangular geometry formed with two base 
reference receivers and one rover station as shown in 
Figure 2 



Fig. 2: Base Reference/Receiver (P) 
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From Figure 2, for a triangle, L = 3, J = 3, h = 0; n = 3(3- 
l)-2x3-0 = 0. This (n = 0) implies the adopting of 
triangular geometry for GPS observations, geometrical 
stability of results is achieved. For a line used as the 
conventional method of GPS relative technique of one 
base reference station and one rover position as shown in 
Figure 3, i.e. computed to be n = -4 to show that line 
used for field observation may not produce stable results. 

A- A 

Base/Ref. L - 1 Rover 

Fig. 3: A line 

For a pentagonal geometry, L = 5, J = 5, h = 0, therefore n 
= 2 to indicate and imply that a deformed and in stable 
shape that may affect the data. 

n. METHOOLOGY 
2.1 GPS Field Observations 

Dual base reference stations were used to determine 
the ellipsoidal heights of the observed controls used as 
rover positions with DGPS receivers and accessories. 
Three online post processing software was used to 
process for the ellipsoidal height and the arithmetic means 
of the ellipsoidal heights were computed. The average 
ellipsoidal heights of each point was used with the 
existing orthometric height to determine the geoid 
undulation of each control point, equation (5). See the 
results in Table 1 


Table 1: Average Ellipsoidal Fleights, Existing Orthometric Fleights and Computed Geoid Undulations 



COORDINATE REGISTER VALUE 

POST 

PROCESSED 
AVERAGE h 

(m) 

UNDULATION, 
N=h-H (m) 

CONTROL 

POINTS 

EASTINGS 
(m) (e)X 

NORTHINGS (m) 
(n)Y 

EXISTING 

ORTHO. 

HEIGHTS 

H (m) 

FCC11S 

331888.114 

998442.043 

485.447 

509.396 

23.949 

FCT260P 

255881.175 

993666.807 

201.944 

224.74 

22.787 

FCT103P 

340639.766 

998375.578 

532.558 

556.836 

24.278 

FCT12P 

333743.992 

1008308.730 

735.707 

760.192 

24.485 

FCT19P 

337452.408 

996344.691 

635.644 

659.824 

24.18 

FCT2168S 

310554.927 

1009739.930 

431.087 

455.274 

24.187 

FCT24P 

322719.776 

1001884.850 

453.804 

477.987 

24.183 

FCT276P 

351983.716 

1025998.314 

625.572 

649.848 

24.276 

FCT4154S 

329953.882 

1003831.280 

476.981 

501.232 

24.251 

FCT4159S 

326124.422 

1003742.860 

452.230 

476.553 

24.323 

FCT66P 

299148.035 

998114.283 

297.111 

321.115 

24.004 

FCT9P 

329821.512 

1007612.091 

497.253 

521.693 

24.440 

FCT35P 

322183.380 

992926.363 

427.171 

451.299 

24.128 

FCT57P 

303234.270 

992916.402 

323.844 

347.795 

23.951 

FCT4028S 

330164.634 

1001388.240 

449.592 

473.942 

24.35 

FCT53P 

308943.361 

993406.773 

351.943 

375.955 

24.012 

FCT4652S 

329441.767 

997474.808 

462.711 

487.113 

24.402 
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FCT162P 

270791.291 

934625.533 

189.696 

215.091 

25.395 

FCT130P 

330982.584 

952889.869 

695.608 

719.383 

23.775 

FCT2327S 

282526.612 

973821.470 

183.287 

207.482 

24.195 

FCT2652S 

271370.273 

945385.429 

138.952 

163.741 

24.789 

FCT2656S 

272644.591 

941062.460 

204.724 

229.229 

24.505 

FCT83P 

332954.205 

987231.606 

568.752 

592.819 

24.067 

XP382 

284074.729 

983364.863 

274.586 

298.390 

23.804 


2.2 Polynomial Surfaces 

The two polynomial surfaces considered to 
represent/model the FCT continuous vertical reference 
surface are: i) multiquadratic and ii) bicubic. 

1. Multi - quadratic model (nine parameters) from 
Sanlioglu et al. (2009) 

N = u n + u.X + u/Y + U-.X ~ + u.Y + 

( 11 ) 

a 5 XY + a 6 X 2 Y + a 7 XY 2 + a s X 2 Y 2 

Multi-quadratic interpolation according to Yanalak and 
Baykal (2001) is an analytical method of representing 
irregular surfaces that involve the summation of quadratic 
surfaces. Kirici and Sisman (2017) stated that even if the 
reference points are not homogeneously distributed, the 
results of surface modelling are barely affected. This is 
particularly applicable to the present studies with 
reference to the lopsided distribution of controls selected 
(after reconnaissance surveys) for use in geometric geoid 
development. 

2. Bi- cubic model (third-order polynomial) 

N = a 00 + qX + Uq^Y + u*,qX ~ + 

u^XY + u Q r,Y + u^qX + ( 12 ) 

a 2l X 2 Y + a n XY 2 +a 0i Y 3 

Where, 

Y = ABS(y- y a ) 

X = ABS (x-x Q ) 

y = Northing coordinate of observed station 
X = Easting coordinate of observed station 
y o = Northing coordinate of the origin (average of 
the northing coordinates) 

X = Easting coordinate of the origin (average of the 
easting coordinates) 

2.3 Least Squares Equation and Solutions 

Observation equation was formed for each point and 
solved to determine the polynomial coefficients X from 
the observation equation generally as given by Ono et al. 
(2004): 

V = AX-L (13) 

where, 


A = coefficient matrix 

X = vector of unknown parameters/coefficients 
L = geoid undulations 

Applying least squares principles, the solution is given by 

X = (A T A) 1 (A r L) for unit weight (14) 

Unit weight (W = 1) is assumed due to equal reliability of 
observations. 

The geoidal undulation of at least six points must be 
known within the study area to enable redundancies for 
the robustness of least squares solution. In this study, 
twenty-four (24) points with both ellipsoidal and 
orthometric heights are known. The model parameters 
determined from least squares solutions are: 

Multi - Quadratic Model Parameters 




' 24.2248901210000000000 ' 



- 0.00002409340580587179 

a 2 


- 0.00008013620770038382 

a 3 


0.00000000009699046795 

a 4 

= 

0.00000000370280953876 

a 5 


0.00000001166702184889 

a 6 


- 0.00000000000021600943 

a 7 


- 0.00000000000045716237 

“sj 


. 0.00000000000000000886 , 


Bicubic Model Parameters 


<0 


' 23.50081592604167925515 N 

a x 


0.00004395285221636346 

a 2 


0.00009105227330487502 

a 3 


-0.00000000156204910634 

a 4 


-0.00000000352164634358 

a 5 


-0.00000000178532065159 

a 6 


0.00000000000004116279 

a 7 


0.00000000000002433215 

a % 


0.00000000000000729517 

a 9 J 


v 0.00000000000002928003 , 


Standard deviation of observations (a) was computed 
using (15): 

\(n- 1) 

where. 
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V = residual = difference between model and known 
geoid height. 
n = number of points. 

2.4 Interpolation of Geoid and Orthometric Height 
Modelling 

Microsoft Excel program was developed to 
interpolate the geoid undulation and hence model the 
orthometric height for each point within the study area. 
The x , y and h are input into the Microsoft Excel program 
developed to interpolate both geoid and orthometric 
heights. The modelled orthometric heights were then 
compared with their corresponding existing orthometric 
heights of the controls and the standard deviation was 
computed from equation (15) as cr lti = 1 I cm and 

^Bicubic = 14cm 


2.5 Hypothesis Testing for Comparison of 


Orthometric Height 

The null hypothesis is given by H 0 

while the 

alternative is H-± and is formulated as follows: 

H„: R = 0, no relationship between H MultiquadEtlc and 

^MSL 

(16) 

H i : R ^ 0, there is a relationship 

between 

^Multiquadatic an d ^MSL 

(17) 


Significance level OC = 0.05 i.e. 95% confidence level. 
Decision rule: reject H a if computed /| > 
from t- distribution table. 


Scenario: H MuMqu!idBtic and H MSL 

The statistical significance of the relationship was 
computed by /-test statistics formula given in 
janda.org/c/10/lectures/topic/L as: 

t = R*J(n-2)/(\-R 2 ) (18) 

In the case of H Multiquadatic and H MS| , the computed t = 
0 while table t =1.717. From the decision rule above, we 
reject H a i.e. the existing heights do not have any 

correlation with the modelled heights. Hence H is 
accepted and further, it may be an indication of 
coincidence of the two surfaces but referenced to different 
reference datum, the geoid and the mean sea level. Height 
values based on the geoid (multi-quadratic or bi-cubic 
models) are the desired orthometric heights and is the 
primary goal of this study in FCT for height 
modernization according to Nwilo (2013). 

2.6 Evaluation of Surface Fitting Techniques 
Alevzakou and Lambrou (2011) stressed the need to 
determine if a surface of higher degree is necessary in 
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geometric geoid modelling by using the relationship 
given as: 

0j(7 2 /r 2 CTi)< F,r 2 (19) 

Where, 

/'i, /' 2 = degrees of freedom of the smaller 
degree surface and the greater surface 
respectively. 

CJj, (J 2 = standard deviation of the two surfaces 
respectively. 

F { ,r 1 =V distribution for one degree difference 
between the tested surfaces CT l = 0.109959231m 
0% = 0.135719119m r x = 15 (multi-quadratic 
surface), r 9 = 14(bi-cubic surface). 

The decision rule is if (/jO% / FCT l ) < F x r~, , then no 
higher surface is needed for geometric geoid modelling of 
FCT. The Computed value of (tjCT 2 / r 2 (J j) = 
1.32242959925 

F ^, r n from F Tables = 4.531 (http://www.stat.ucla.edu/- 
dinov). 

From the relationship (/j<T 2 / 1' 2 (J ] ) and F distribution 

7^,r 0 , since 1.32242959925<4.531, no higher degree 
surface is needed for geometric geoid model in the FCT. 
This is an indication that either multi-quadratic or bi¬ 
cubic model can be used to model orthometric height 
although the multi-quadratic model performed better and 
could be taken as the optimum. 

2.7 Skill and Bias Estimates 

The skill parameter can be used as a measure of the 
model predictive capacity in relation to the observations. 
This skill parameter ranges from negative values to one 
with the corresponding value of one implying a total 
agreement between observations and the model results. 
The bias values computed as zero simply imply that the 
data used, equipment used and personal error did not 
show any bias whatsoever in this study. Bias and skill 
were calculated by the relationship given by Sutherland 
et al. (2004) 

This also suggested that the selection and 
combination of equipment, personnel, field techniques 
and processing methods adopted yielded high quality data 
to produce the FCT geoid surface information as much as 
possible. Orthometric heights from the surface are hence 
based on geoid and compatible with GNSS technique and 
the adopted dual base reference stations technique. 

2.8 Diagnostic Test for Multiquadratic and Bicubic 
Models 
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To carry out a diagnostic test for the predictive 
ability of the models in orthometric height modelling as 
stated in Sinha and Prasad (1979), the computation was 

carried out using 1.98 /J~N where N is the number of 
stations = 24 in this study to compare with chi-squares 
table values. The decision rule is if 1.98 /y[~N < , 

then models are satisfactory at 95% confidence level. In 
this study, 1.98/V/V =0.404. Using the Chi squares 
(% 2 ) test at the 95% degrees of freedom (d.o.f), we have 

for multiquadratic model, degree of freedom = 15, at 

2 

95%, = 24.996; bicubic model, degree of freedom = 

14 at 95% = 23.685. Since 0.404 < 24.996 or 23.685, 

the models proved satisfactory at 95% confidence limits 
for modelling orthometric heights from GNSS techniques 
as confirmed by the diagnostic tests. 

2.9 Application Areas and Importance of the Geoid 
Model 

Applications of geoid are: 

1. For transforming GPS ellipsoidal heights (h) to 
orthometric height for practical surveying and 
engineering applications. 

2. An important part of a National Geodetic Data 
Infrastructure (NGDI). 
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3. The geoidal map can also be used to interpolate for 
geoid heights at any point of interest in 
FCT. 

4. This is useful where the conventional method of 
spirit levelling is costly, tedious, time-consuming 
and costly especially in highly urbanized areas. 

The importance of the determined geoid models in 
orthometric height derivation are:- 

1. Consistency and compatibility with GNSS technique 
is achieved with these models for orthometric height 
determination. 

2. Orthometric heights can be interpolated for all 
points of interest within the FCT. 

III. DISCUSSION OF RESULTS 
3.1 Plot of Geoid Undulation Against Controls 

Figure 5 presents the plots of multi-quadratic and bicubic 
geoid heights against control stations. This was done to 
show graphically the differences between the multi¬ 
quadratic and the bi-cubic models' geoid undulations. It 
can be seen that the two surfaces are nearly coincident 
and identical from a visual inspection of Figure 5. This 
implies and confirms the interchangeability and 
acceptability of the two models for orthometric height 
determination within the FCT but attaching more weight 
to the multi-quadratic model. However, visual 
observations/interpretations are generally subjective. 



Fig. 5: Geoid Undulation Models from both Multi-Quadratic and Bi-Cubic Models. 


3.2 Similarity of Surfaces 

Figure 6 shows the plots of the multi-quadratic model, 
bicubic model and existing orthometric heights. This was 
also done to show graphically the similarity of the three 


surfaces. From Figure 6, it can be confirmed that the 
multi-quadratic model is more suitable and adequate as it 
is smoother for the modelling of orthometric heights in 
FCT by GNSS technique. 
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Fig. 6: Similarity Plots of the Three Surfaces 


IV. CONCLUSIONS 

From the study, the following conclusions can be made: 

1. Multiquadratic model takes good care of the lack of 
homogeneous distribution of selected controls in 
geoid modelling. 

2. Coefficient of correlation (R) and coefficient of 

determination ( R) values of R = 0.995m and 

(R 2 ) = 99% respectively indicate the multi¬ 
quadratic model has a high predictive ability at 95% 
confidence limits. 

3. Dual base reference stations were adopted for data 
acquisition instead of conventional single base 
reference station. This enabled the stability of results 
by exploiting Grashof’s law of stability of triangles. 

4. The feasibility of developing a geoid model for GPS 
user community by GNS S/Levelling in FCT has 
been demonstrated as an alternative approach to 
conventional spirit levelling in orthometric height 
determination 
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